(N 

o 

(N 



43 



1/3 

ll; 

1/3 ' 

o 

</3 . 
>v 
43 ' 
Of 



> 

cn 

in 
o 

(N 



X 
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A multi-chord fiber-coupled interferometer [Merritt et al., Rev. Sci. Instrum. 83, 033506 (2012)] is being used 
to make time- resolved density measurements of supersonic argon plasma jets on the Plasma Liner Experiment 
[Hsu et al., Bull. Amer. Phys. Soc. 56, 307 (2011)]. The long coherence length of the laser (> 10 m) allows 
signal and reference path lengths to be mismatched by many meters without signal degradation, making for a 
greatly simplified optical layout. Measured interferometry phase shifts are consistent with a partially ionized 
plasma in which an initially positive phase shift becomes negative when the ionization fraction drops below a 
certain threshold. In this case, both free electrons and bound electrons in ions and neutral atoms contribute 
to the index of refraction. This paper illustrates how the interferometry data, aided by numerical modeling, 
are used to derive total jet density, jet propagation velocity (~ 15-50 km/s), jet length (~ 20-100 cm), and 
3D expansion. 



I. INTRODUCTION 

Supersonic plasma jets are being studied on the Plasma 
Liner Experiment (PLX)~ as a method of forming im- 
ploding spherical plasma liners^ which are a potential 
standoff driver for magneto-inertial fusion.— The plasma 
jets are formed and accelerated by railguns developed and 
built by Hyper V Technologies and have initial densities 
on the order of 10 17 cm -3 , radii of a few centimeters, 
lengths of a few tens of centimeters, and velocities up 
to 50 km/s4 An eight-chord, visible, heterodyne, fiber- 
coupled interferometer— has been designed and built for 
measuring the plasma density of a propagating single jet 
and the merging of multiple jets. In single jet propa- 
gation studies thus far, interferometry phase shift data 
show a positive value (due to free electrons) becoming 
negative (due to bound electrons associated with neu- 
trals and ions) later in time, consistent with the ioniza- 
tion fraction / dropping as the jet propagates and under- 
goes adiabatic expansion and radiative cooling. In such 
partially ionized plasmas, a two-color interferometer^ 
could simultaneously constrain / and the ion plus neutral 
density n to t- However, our single color system requires 
independent information on either / or n tot to constrain 
the other. 

The main purpose of the paper is to describe the inter- 
pretation of interferometer phase shift signals from par- 
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tially ionized argon plasmas, to show the methodology 
of deducing plasma jet properties from the phase shift 
data, and to illustrate how synthetic interferometer data 
from plasma jet simulations can be used to help inter- 
pret and constrain the experimentally-deduced jet prop- 
erties. More detailed physics results on the plasma jets 
relating to the PLX project will be reported elsewhere. 
The remainder of the paper is organized as follows. Sec- 
tion II very briefly describes the interferometer system. 
Section III derives the interferometer phase shift formula 
as a function of / and ntot for argon plasmas. Section IV 
illustrates how several jet physical quantities are deduced 
from the phase shift data, and Sec. V shows synthetic 
interferometer phase shift data that help constrain the 
interpretation of the experimental data. Section VI pro- 
vides a summary. 



II. INTERFEROMETER DESIGN AND SETUP 

The details of the eight-chord interferometer system 
are covered in a previous publication^ The interferom- 
eter uses a 320 mW, Ao = 561 nm laser with a 10 m 
coherence length, as specified by the manufacturer. The 
long coherence length of the laser allows for a sub-fringe 
phase resolution of 4° even with a path length mismatch 
between the probe and reference beams of a few tens of 
centimeters arising from the setup geometry. Figure 1 
shows the experimental setup, in which the eight probe 
beams are positioned approximately perpendicularly to 
the direction of jet propagation (Z axis), starting at Z = 
35 cm in 6.35 cm intervals. 



2 




FIG. 1. Schematic of the experiment in a 9 ft. diameter 
chamber. Eight interferometer chords each separated by 6.35 
cm intersect the path of plasma jet propagation (Z axis) as 
shown. The plasma jet is launched at Z—0, and the first chord 
is at Z=35 cm. 



III. PHASE SHIFT ANALYSIS FOR A PARTIALLY 
IONIZED PLASMA 



The index of refraction is related to the phase shift mea- 
sured by the interferometer as follows: 8 



A4>=^f{N-l)dl = ^f 



SNdl, 



(4) 



where A is the wavelength of the laser. The phase change 
due to neutral atoms in the plasma is then^ 
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Like the index of refraction of a neutral gas, the index 
of refraction for ions is also a result of the perturbation 
of the incident light by bound electrons in the ion. So, 
the Dale- Gladstone relation is still valid for the ions in 
the plasma. The phase change due to ions in the plasma 
is 



2tt 

A(j>i = —l\,m, I ii} ill. 



(6) 



For a strongly ionized plasma, the index of refraction 
for a laser beam of wavelength Ao can be approximated 
as a function of only the line-integrated electron density 
of the plasma: 8 



e 2 A 



47reorn e c 2 



n e dl. 



(1) 



In this approximation it is assumed that the density of 
neutral particles in the plasma is negligible and that the 
frequency of the laser light has been chosen such that any 
interaction with the ions is also negligible. However, in 
a partially or weakly ionized plasma these assumptions 
are no longer valid and the ion and neutral particle den- 
sities may both contribute significantly to the index of 
refraction of the plasma. 



A. Neutral Atom and Ion Contributions to the Phase 
Change 

The index of refraction of a gas can be found using the 
Dale-Gladstone relation^ 



Kp = (N - 1) = 5N, 



(2) 



where p is the mass density of the gas, K is the specific re- 
fractivity, N is the refractive index, and SN = N— 1. The 
specific refractivity of the gas is determined by the gas 
species and ionization state. For most neutral gases, the 
index of refraction and the density are well documented 
at standard temperature and pressure (STP) and can be 
used to calculate the specific refractivity. Thus the index 
of refraction of the neutral atoms in the plasma is&i^ 



5N n = Km n n 7 , 
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B. Expression for the Total Phase Shift 

The total interferometer phase shift is a superposition 
of the electron, ion, and neutral particle contributions 
since the indices of refraction are additive^ The interfer- 
ometer signal analysis has been defined such that the 
contribution from any material with N < 1 will give 
Acj) > 0, and materials with N > 1 will give Acf> < 0. 
Since N e < 1 and Ni,N n > 1, then the total phase shift 
is : 



Acj> tot = A0 e - A0 
Ae 2 
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Because the phase change from the ions and neutrals is 
in the opposite direction of the phase change from the 
electrons, using the approximation A<j)tot — A<j) e can un- 
derestimate the actual electron density of the plasma. 

The plasma jets are cold (T e w leV) by the time they 
intersect the interferometer chords, and thus the plasma 
is assumed to consist of only singly-ionized and neutral 
atoms, i.e., n e = ni. The total phase change is thus 



A<t> t . 
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which is a function of only the line-integrated ion and 
neutral densities. 
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Phase Shift as a Function of Ionization Fraction 
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FIG. 2. Calculate phase shift versus ionization fraction for an 
argon plasma at a variety of PLX-relevant atomic densities. 
The phase shift of A<j> to t = occurs at f = 0.072350. The 
gray areas are inaccessible parameters. 



The phase shift from singly ionized argon is then 
P Ae 2 



47reom e c 2 



/ 



-(1 + 0.67/-/)^^^-] / n tot dZ (13) 
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= 1.6204 x 1CT 17 [/ - 0.07235] / ra tot dZ, (14) 

a function of only / and n to t- For argon, A<f> to t = 
occurs for / = 0.07235. The most general form for the 
phase shift is then 



Afoot = C[f -f ] / n tot dl, 



(15) 



where C is a scale factor for the gas species and fo is the 
ionization fraction for which Afoot = 0. Figure 1 shows 
a plot of A(f> t ot versus / for a variety of J ntotdl values 
of relevance for a single plasma jet. Since J ntotdl is 
always positive, then A<f> > requires f > fo and A<j> < 
requires f < f . 



C. Calculating the Total Atomic Density 



IV. EXPERIMENTALLY DERIVABLE PARAMETERS 



The total atomic density of the plasma is defined as 
the sum of the ion and neutral densities, n to t = + n n . 
The interferometer is not equally sensitive to ion and 
neutral densities so a second relation between rij and n n 
is required in order to determine the total phase shift in 
terms of total atomic density. Assuming a value for the 
ionization fraction, 



/ = 



71, 



n t + n n n to t ' 



(9) 



fulfills this requirement. The total phase shift can now 
be written in terms of / and the total atomic density ntot- 
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Assuming an approximately uniform ionization along the 
path of the laser through the plasma, then the phase shift 
simplifies to 

A(ptot = [{-. o ~ -^ K i m i)f 

-^^S^(l - /)] / ntotdl. (11) 
A n„ J 



D. Parameters with Argon 



PLX is currently studying argon single jet propagation 
from a HyperV Mark I plasma railgun 4 with < 350 kA of 
peak current. Figure 3(a) shows the phase shift (left axis) 
versus time averaged over 25 shots with I pea k = 260 ± 5 
kA and a resulting chamber pressure rise after each shot 
of P = 0.4±0.1 mtorr. The signal for the averaged 
jet displays a sharp rise at the beginning of the signal, 
multiple peaks, and a transition from Acf> > to A<p < 0. 



A. Ionization Information 

Since the phase shift is a function of both / and n tot , 
then independent information, such as from spectroscopy 
or synthetic diagnostics, is required to estimate / and 
determine J ntotdl (or vice versa). However, the sign of 
A<j>tot alone can bound /. Any signal with A(f> t ot > 
has a lower bound on the ionization fraction of fo. As- 
suming a uniform / along the path of the laser through 
the plasma, then according to Eq. (15) any signal with 
Aifitot < has f < fo- In Fig. 3(a) the peak at 20 fxs < t < 
60 us must have / > 0.07235 and all subsequent peak for 
60 fis < t < 150 (as must have / < 0.07235. Thus the ion- 
ization is not constant along the jet's axial length. The 
ionization bounds set by the experimental signal can be 
used to calibrate simulated jet initial conditions such as 
initial jet temperatures, densities, and temperature and 
density distributions. 



For argon, the Slater screening constant gives a specific 
refractivity for singly ionized argon>ii 

Katii = 0-67 x KatI = 0-67 x — — . (12) 

niArnArl 



B. Total Atomic Density and Jet Structure 

Once independent information about the ionization 
fraction ia available, then the interferometer can deter- 
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FIG. 3. (a) Phase shift versus time for the eight chords (distances from plasma gun indicated in legend). The /-dependent 
line-integrated density is given on the right-hand y-axis. Arrival times of the leading edge and trailing peaks are marked as 
vertical lines, (b) Line-integrated ion plus neutral density versus time assuming f ~ fo + 0.01 for A(j> > and f = f — 0.01 
for A<f> < 0. Times of the leading and falling edges (half peak) are shown by vertical lines, (c) Peak line-integrated density 
versus time for the peaks indicated in (b). (d) Jet length versus time for the leading (left y-axis) and trailing jet (right y-axis) 
structures. 



mine the total atomic density as well as several related 
jet parameters. 

For purely illustrative purposes we assume an arbitrary 
ionization fraction of the form f(t) = fo + a for A(f> > 
and f(t) = fo — a for A<fr < 0, where a = 0.01. This func- 
tional form was chosen because it is the simplest that ful- 
fills the requirements that f > fo for positive phase shifts 
and / < fo for negative phase shifts. The actual f(t) for 
the jet is currently under investigation. Figure 3(b) is a 
plot of J ritotdl versus time for the averaged jet in Fig. 
3(a) assuming the arbitrary form of f(t). The first dis- 
tinguishable parameter is the order of magnitude of the 
line-integrated jet atomic density, / n to tdl ~ 10 18 cm~ 2 
in Fig. 3(b). Higher jet densities are desirable for creat- 
ing plasma liners that reach high pressure^ because the 
liner ram pressure scales linearly with liner density. 

These jets have a consistent set of distinguishable 
structures: a sharp rise in atomic density that marks the 



leading edge of the jet, two distinct jet sections separated 
by a region of very low density, a single peak in the first 
jet section, and a double peak in the second jet section. 
The multiple peaks in the jet are possibly due to ringing 
of the gun electrical current. 



C. Jet Velocity 

Distinct structures in the signal have distinguishable 
arrival times at each interferometer chord. The distance 
between each chord is 6.35 cm, so the velocity of each dis- 
tinguishable structure can be calculated using the time 
difference between the structure arrival times for consec- 
utive interferometer chords. The jet velocity at the front 
of the jet is calculated using an arrival time of the lead- 
ing edge calculated at half the peak value of the leading 
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peak, t = t(A$ e ak /2) as shown in Figs. 3(a)-(b). The 
leading edge of the jet has a velocity of v^ 9 « 54 km/s 
(averaged over 25 shots). The velocity for later peaks 
in the jet can also be calculated; the peak velocity for 
the final jet peak is calculated using the time difference 
between the peaks themselves for consecutive interferom- 
eter chords. The final peak velocity is v"^ 9 ~ 15 km/s. 



D. Jet Expansion 

Phase shift data can be used to estimate jet radial 
and axial expansion as a function of time by invoking 
conservation of mass: 

Utotif) x V(t) — constant 
(n tot D)(t) x D(t) x L(t) — constant, (16) 

where V(t) is the volume of the jet, 
D(t) 



(n t otD)(t) x L(t) 



(17) 



is the jet diameter, L(t) is the jet length, and J n tot dL 
(ntotD)(t). Jet length is determined by 



L 



jet 



J 



-(tf-ti). 



(18) 



where i and / denote the initial (leading) and final 
(falling) edges of each sub-jet, respectively, and the ar- 
rival times of the edges are taken at half the jet peak 
atomic density as indicated in Fig. 3(b). Jet length ver- 
sus time for both sub-jets is plotted in Fig. 3(d). L(t) is 
non-constant for both sub-jets and implies axial expan- 
sion of the jet. The value of (n to tD)(t) for a sub-jet is 
determined from the maximum value and arrival time of 
a peak of the sub-jet. Fig. 3(c) plots (/ n tot dl) peak (t) vs 
time for the peak indicated in Fig. 3(b) with quadratic 
and exponential fits to the data. These equations can be 
used to determine 



ntot(t) = 



(n t otD)(t) 
D(t) 



oc (n tot D) 2 {t) x L(t). (19) 



in terms of experimentally-derived parameters. 



V. SYNTHETIC INTERFEROMETRY DATA 

For synthetic interfcromctry, the phase shift formula, 
Eq.(14), has been applied to simulated jet data from the 
single-fluid Nautilus^ and two-fluid Large Plasma Sim- 
ulation (LSP) 1 - codes. Chord locations for the synthetic 
interferometry match those in the experiment. Relative 
positive and negative phase shift amplitudes, jet struc- 
tures such as the sharp leading jet edge and the transi- 
tion between positive and negative phase shifts, and jet 
arrival times at the interferometer chords are all used to 
calibrate simulation conditions. 



Figures 4(a) and 4(b) show synthetic interferome- 
ter phase shift signals from ID LSP— two-fluid sim- 
ulations and 3D Nautilus^ single-fluid simulations, 
respectively. Both simulations include non-local- 
thermodynamic-equilibrium (non-LTE) equation-of-state 
(EOS) tabular models^ and assume a leading and a 
trailing jet with initial conditions as given in Fig. 4. 
The ID LSP simulation assumes that the radius of both 
jets expands thermally for its path length calculation. 
The trailing jet requires a lower temperature and cor- 
responding / to generate a negative instead of positive 
phase shift. For these conditions both the negative and 
positive signals have roughly equal amplitudes, and are 
in approximate agreement with the experimental signal 
(Fig. 2). The non-LTE EOS model used bounds the jet 
temperature at T e < 1 eV for / <C 1. The leading edges of 
the synthetic phase shift data also show a fast rise con- 
sistent with the experimental data. These simulations, 
which have produced synthetic data in good agreement 
with the experimental data, will be used to interpret and 
better understand the physics of single jet propagation. 



VI. SUMMARY 

An eight-chord, fiber-coupled interferometer with a 
long coherence length laser— is being used to study su- 
personic argon plasma jets. The long coherence length 
allows mismatches between probe and reference beams, 
allowing the use of only one reference beam for all eight 
probe beams, and simpler reconfigur ability of the chords 
for doing different experiments. In single jet propaga- 
tion experiments, the interferometer phase shift shows 
a positive phase shift becoming negative later in time. 
This is shown to be consistent with the / of the plasma 
jet falling while the jet propagates. We derive a general 
phase shift formula accounting for free electrons, neutral 
atoms, and singly ionized argon ions. The formula de- 
pends only on n to t and /. Examples are shown for how 
the interferometer data, with some assumptions or inde- 
pendent input from other diagnostics or numerical mod- 
eling, can be used to deduce jet physical properties such 
as ritoti /) velocity, and jet expansion. Finally, synthetic 
interferometer phase shift data calculated using Eq. (14) 
from single- and two-fluid simulations of argon plasma 
jet propagation are shown, giving an example of how the 
synthetic data can help with interpretation of the ex- 
perimental phase shift data. Physics conclusions of the 
plasma jet studies as they relate to the PLX project will 
be reported elsewhere. 
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FIG. 4. Synthetic interferometer phase shift data from (left) LSP two-fluid and (right) Nautilus single-fluid jet simulations. 
Initial conditions used in the simulations are indicated in the plots. Synthetic chord positions relative to the plasma gun are 
given in the legends. 
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